
/* 
robustness to additional controls
*/

set more off

estimates clear
clear

use ".\Dropbox\Vaccination\covid.dta" // load the data

drop if statefips==48|statefips==15 // drop states without the vaccination data

zscore(rep_pres2016 vacc*) //normalizing the variables

label variable z_rep_pres2016 "Trump vote share, 2016"


// Trump vote share

reg z_vaccinated_18pluspop_ratio z_rep_pres2016 conservative lat lon  temp rain income native_share white_share col_share elder_share male_share  mfg_share_diff popden i.statefips, cluster(statefips) 
estimates store trump_conservative

reg z_vaccinated_18pluspop_ratio z_rep_pres2016 republican conservative lat lon  temp rain income native_share white_share col_share elder_share male_share  mfg_share_diff popde i.statefips, cluster(statefips) 
estimates store trump_republican

reg z_vaccinated_18pluspop_ratio z_rep_pres2016 tox conservative lat lon  temp rain income native_share white_share col_share elder_share male_share  mfg_share_diff popde i.statefips, cluster(statefips) 
estimates store trump_svi

reg z_vaccinated_18pluspop_ratio z_rep_pres2016 healthenvironmenthealthdispariti  conservative lat lon  temp rain income native_share white_share col_share elder_share male_share  mfg_share_diff popde i.statefips, cluster(statefips) 
estimates store trump_disparity

reg z_vaccinated_18pluspop_ratio z_rep_pres2016 sk2014 conservative lat lon  temp rain income native_share white_share col_share elder_share male_share  mfg_share_diff popde i.statefips, cluster(statefips) 
estimates store trump_sk

reg z_vaccinated_18pluspop_ratio z_rep_pres2016 flu conservative lat lon  temp rain income native_share white_share col_share elder_share male_share  mfg_share_diff popde i.statefips, cluster(statefips) 
estimates store trump_flu

reg z_vaccinated_18pluspop_ratio z_rep_pres2016 conservative republican tox flu healthenvironmenthealthdispariti  sk2014  lat lon  temp rain income native_share white_share col_share elder_share male_share  mfg_share_diff popde i.statefips, cluster(statefips) 
estimates store trump_all

estout trump* using ".\Dropbox\Vaccination\Draft\tab_robust_ols.tex", replace style(tex) cells(b(fmt(3) star) se(par fmt(3))) stats(r2 N, fmt(2 0) labels("R-squared" "Observations")) keep(z_rep_pres2016) label mlabels(none) collabels(none) starlevels(* 0.10 ** 0.05 *** 0.01)
